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Abstract: 

Gamma-emitting radionuclide ”"Tc is globally used for the diagnosis of various pathological 
conditions owing to its ideal single-photon emission computed tomography (SPECT) 
characteristics. However, the short half-life of °"Tc (Ti2 = 6h) makes it difficult to store or 
transport. Thus, the production of ”Tc is tied to its parent radionuclide ”Mo (Ti. = 66h). The major 
production paths are based on accelerators and research reactors. The reactor process presents the 
potential for nuclear proliferation owing to its use of highly enriched U (HEU). Accelerator-based 
methods tend to use deuterium-trittum (D-T) neutron sources but are hindered by the high cost of 
tritium and its challenging operation. In this study, a new Mo production design was developed 
based on a deuterium-deuterium (D-D) gas dynamic trap fusion neutron source (GDT-FNS) and a 
subcritical blanket system (SBS) assembly with a low-enriched U (LEU) solution. GDT-FNS can 
provide a relatively high neutron intensity, which is one of the advantages of Mo production. We 
provide a Monte Carlo-based neutronics analysis covering the calculation of the subcritical 
multiplication factor (k,) of the SBS, optimization design for the reflector, shielding layer, and *°Mo 
production capacity. Other calculations, including the neutron flux and nuclear heating distributions, 
are also provided for an overall evaluation of the production system. The results demonstrated that 
the SBS meets the nuclear critical safety design requirement (k,<0.97) and maintained a high °Mo 
production capacity. The proposed system can generate approximately 157 Ci Mo for a stable 24 h 
operation with a neutron intensity of 1<10'* n/s, which can meet 50% of China’s demand in 2025. 


Key words: Gas dynamic trap, Fusion neutron source, Molybdenum-99, Low-enriched uranium, Subcritical 
blanket system 


1 - Introduction 


According to the World Nuclear Energy Association (WNA) [1], over 10000 hospitals 
worldwide use radioisotopes for the diagnosis and treatment of diseases. It is estimated that 
approximately 90% of radioisotopes are used in diagnostic procedures, among which Tc is the 
most commonly used, which benefits from its ideal characteristics of single-photon emission 
computed tomography (SPECT). Further analysis also demonstrates that the use of *"Tc constitutes 
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approximately 80% of nuclear medicine procedures, while 85% are used for diagnostic scans 
(Updated April 2022) [1]. However, the production of ®™Te (Ti2 = 6 h) is tied to its mother 
radionuclide Mo (Ti = 66 h), which makes the study of Mo production an essential research 
topic. 

The production of ”"Tc consists of two steps: 1) production of Mo via mechanisms indicated 
in the schematic diagram of Fig. 1 (neutron-fission, gamma-fission, neutron-gamma, gamma- 
neutron, etc.) and 2) separation of °™Tc after it decays from Mo through the emission of beta 
particles [2]. As indicated, °Mo and *™Tc have short half-lives. Nonetheless, unlike the relatively 
short half-life of °"Tc (Tı2= 6 h), which is the main hindrance for its transport, the half-life of °Mo 
(Ti2= 66 h) allows a relatively adequate time for transportation. Based on the half-life and transport 
time, the kinetics of producing "Te depend on the production of *Mo; that is, under normal 
circumstances, the resulting product Mo is transported to the target country immediately after 
production, and subsequently its decayed product (*"Tc) is quickly extracted and used in hospitals 
or nuclear medical centers. 
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Fig. 1 Schematic diagram of the °™Tc production process 


The global shortage of Mo has recently been mainly attributed to aging nuclear reactors and 
their decommissioning [3-5]. A typical example is the event of the National Research Universal 
(NRU) in Canada, which produced approximately 40% of the world’s Mo supply but ceased 
production by October 31“, 2016 [3]. Currently, most Mo isotopes are produced by the fission 
process with highly enriched U (HEU) reactors. These nuclear reactors are the Belgian Reactor 2 
(BR-2), High Flux Reactor (HFR) in the Netherlands, LVR-15 Reactor in the Czech Republic, 
Maria Research Reactor (Maria) in Poland, Open Pool Australian Light Water Reactor (OPAL), and 
South African Fundamental Atomic Research Installation (Safari-1) [6]. Additional information 
regarding *’Mo production, specific reactors, and target materials are listed in Table 1. Most of these 
reactors are not only facing the problems of aging and decommissioning, but also pose great risks of 
nuclear proliferation. 


Table 1 Main *Mo production reactors[6] 


Country Reactor Targets Capacity Start Est. stop 
Belgium BR-2 HEU 6500 1961 2026 
Netherlands HFR LEU™ 6200 1961 2022 
Australia OPAL LEU 3500 2006 2030 
Czech Republic LVR-15 HEU 3000 1989 2028 
South Africa Safari-1 LEU 3000 1965 2025 
Poland Maria LEU 2200 1974 2030 
Russia RIAR: three HEU 890 1961-70 


USA MURR HEU 750 1966 
Argentina RA-3 LEU 400 1967 2027 
Total 26,440 


* 6-day Ci ”Mo: a 6-day Ci is the measurement of the remaining radioactivity of ”Mo six days after it leaves the processing 
facility (end of processing (EOP)) 
* In March 2021 HFR started using LEU for the production of radioisotopes. 


To solve the problems of ”Mo shortage and the nuclear proliferation of HEU targets, scientists 
have proposed a series of new methods to replace the path of HEU fission to produce Mo. These 
methods can be divided into three categories: 1) **°U(n, f)”Mo reaction in low-enriched U (LEU) 
reactors [7, 8]; 2) solid target irradiation based on an accelerator, such as the neutron capture 
Moln, y)’’Mo reaction [9-13], '°’Mo(n, 2n)”’Mo reaction [14,15], Mop, 2n)”Mo reaction [16], 
'Mo(y, n)’”’Mo reaction [17], and photon-induced reaction of **U fission **U(y, f)”Mo [18]; and 3) 
LEU solution fission of the **°U(n, f)’?Mo reaction in subcritical systems [19-23]. Among these, the 
last method is the most efficient and reliable means of production and has become a prime choice 
for ”Mo production owing to the following advantages:1) compared to the solid target irradiation 
method based on an accelerator, the LEU solution can be effectively recovered and reused, thus 
significantly reducing the generation of radioactive waste; 2) the LEU fission method has a high 
production efficiency and low cost; and 3) comparing the LEU solution in subcortical systems to 
the HEU fission method, there is an advantage of avoiding nuclear critical safety accidents, and it 
also prevents nuclear proliferation. It is also relatively easy to apply for a license for construction 
and operation. 

In 2021, Han et al. [22, 23] proposed a subcritical **Mo production system driven by an 
accelerator-based D-T neutron source, in which the accelerated deuterium ions bombard the tritium 
target, and the deuterium-tritium (D-T) fusion reaction generates neutrons. The LEU solution target 
is then irradiated by neutrons for the fission of **°U via the **U(n, f)”Mo reaction. Although this 
method does not require a supply of HEU, it has the disadvantage of using tritium. In addition to its 
high cost, it is difficult to obtain licenses for owning and operating tritium. 

In this study, we propose a new design for ”Mo production based on the LEU subcritical 
blanket system (SBS). The system is driven by a gas dynamic trap-based fusion neutron source 
(GDT-FNS). However, instead of the normal deuterium-tritium fusion reaction, a deuterium- 
deuterium (D-D) fusion reaction neutron is used to induce fission in *°U. In addition to the apparent 
advantage of avoiding the HEU system, our proposed *Mo production system has numerous 
advantages, such as a compact structure, high neutron source intensity, and lack of tritium 
consumption, leading to low capital costs. To ensure the design safety and process optimization for 
the SBS ”Mo production path, a neutronics analysis of the production system was performed using 
the SuperMC code (Monte Carlo Particle Transport code). A detailed neutronics analysis is provided 
for an overall evaluation of the production system, including the analysis of the subcritical 
multiplication factor (k,), neutron flux, and heat deposition. A set of optimization parameters of the 
production system was obtained by maintaining a relatively high production rate and safety 
standard, including the geometric size, material components, and concentration of the LEU solution. 


2 Model and Method 


The LEU solution SBS driven by GDT-FNS mainly includes GDT-FNS and a *Mo SBS, as 


shown in the schematic diagram in Fig. 2. The SBS is arranged in the high neutron flux region of 
the GDT, forming a fan-shaped blanket structure. Detailed descriptions are provided in sections 2.1 
and 2.2. 
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Fig. 2 Schematic diagram of the LEU solution SBS driven by GDT-FNS 


2.1 Gas Dynamic Trap-based Fusion Neutron Source (GDT-FNS) 


A GDT is a type of axisymmetric magnetic mirror device [24, 25]. Under the action of a 
specific magnetic field, the warm plasma constrained in the GDT vacuum chamber frequently 
collides and causes a fusion reaction, which can provide D-D or D-T fusion neutron sources, and 
the neutrons are "high at both ends and low in the middle.” This type of neutron source has the 
advantages of a high neutron flux, large testing space, compact structure, and low construction cost. 
The GDT-FNS not only meets the performance of fusion materials/components, but its resulting 
high neutron flux can be used to conduct the study of applied nuclear technology, such as medical 
isotope production, neutron photography, neutron irradiation (breeding), and for the low-dose 
neutron effect in cells. 

In this study, the GDT-FNS designed by the Hefei Institute of Physical Science [26, 27] was 
used to analyze the neutronics of a solution-based-LEU SBS. Quasi-monoenergetic neutrons with 
an energy of approximately 2.5 MeV were generated by means of the D-D () reaction in the central 
vacuum system (CVS). As indicated in Fig. 2, a GDT-FNS mainly includes a neutral beam injection 
system (NBI), CVS, magnetic coil (MC), and a high neutron flux region. The main parameters of 
the GDT-FNS with the D-D operation are listed in Table 2. Owing to the axisymmetric 
characteristics of GDT-FNS, SBSs with flexible settings can be arranged in the high neutron flux 
region to meet the increasing demand of ”Mo production. 


Table 2 Main parameters of GDT-FNS with D-D operation 


Parameters Value 
Total length of GDT-FNS 2000 cm 
Length of high neutron flux region 100 cm 
Plasma tube radius at high neutron flux region 35 cm 


Vacuum vessel radius at high neutron flux 100 cm 


region 


Magnetic field in center 0.76 (0.94") T 


Magnetic field in plug 25.8 T 
NBI power 50 MW 
NBI angle 45 degree 
NBI energy 70 keV 
Gas feed rate 4200 eq.A 
Peak fast ion density 7x10” m” 
Bulk plasma density 5x10” m” 
Electron temperature 1150 eV 
Bulk ion temperature 2340 eV 
Total neutron intensity of GDT-FNS 9.04x10" n/s 


Neutron intensity of high neutron flux region 8.54x10'*n/s 
* At the start of discharge when the plasma beta is near zero. 


The plasma parameters of the GDT-FNS system were simulated using the 1-D code DOL [28], 
which is based on a nonstationary numerical model describing the confinement of two different 
plasma components. During the simulations, a pure deuterium beam was injected into the central 
vacuum chamber, and the axial distributions of the D-D neutron generation rate were obtained, as 
shown in Fig. 3. The results demonstrate two high neutron flux regions, which are between the 
-700—600 cm and 600-700 cm axial positions of the GDT-FNS. This important finding is the 
main reason for arranging the SBS in the high neutron flux region. 


Fig. 3 The axial distribution of the D-D neutron generation rate of GDT-FNS 


2.2 Subcritical Blanket System (SBS) 


The SBS for the production of Mo is arranged in the high neutron flux region of the GDT- 
FNS and forms a fan-shaped blanket structure. The main material components of the SBS include 
the LEU solution, solution container, reflector, and shielding layer. A schematic of the SBS Mo 
production model is shown in Fig. 4. Based on the preliminary analysis and considering the FNS 
design constraints, certain geometrical and material parameters were fixed, such as the dimensions 
and specific material composition. The SBS was 100 cm long and less than 100 cm thick. The 
thickness of the LEU solution ranged between 30-50 cm, and the solid angle of the LEU solution 
relative to the central axis of the neutron source was between 2/4-2/3. Thus, the variable parameters 
that need to be optimized for this study include the thicknesses of the LEU solution, reflector, and 
shielding, and the material types for the reflector and shielding layer. Considering the LEU solution, 
the UO.SO;, solution was selected with varying U concentrations (60 g/L to 150 g/L), with a *°U 
enrichment of 19.75%. The low margin (60 g/L) was selected to obtain a considerable output of 
Mo, and the upper limit (150 g/L) was based on the saturated U (UO2SO,) concentration at room 
temperature. 
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Fig. 4 Schematic structure diagram of SBS 


2.4 Calculation Method 


2.4.1 Calculation program and uncertainty 


The neutronics parameters of the SBS were calculated using the Super Monte Carlo Simulation 
Program (SuperMC) version 3.2 [29] coupled with ENDF-VII cross-section libraries. Herein, the 
steady-state neutronics parameters of the SBS mainly included the ks, neutron energy spectrum, 
activity of the produced ”Mo, and heat deposition. In this study, 10 million statistical particles were 
used for each neutronics calculation, and the corresponding statistical uncertainties were less than 
1%, except for the calculations of the energy deposition, where the statistical uncertainties were less 
than 3%. 


2.4.2 Subcritical multiplication factor 


The critical safety state of a subcritical system can be characterized by the k, when the U 
fission system has an externally driven neutron source. Parameter k, is defined as the ratio of the 
fission neutron number to the total neutron number in the system [30,31], as shown in equation (1). 

(1) 
Here, is the intensity of the externally driven neutron source [n/s], is the fission reaction rate 
[fission/s], and is the average number of neutrons generated by a fission reaction. The parameter ks 
is usually required to be less than 0.98 to ensure the operational safety of a subcritical system [32- 
36]. The calculation of ks was performed with a general source card (SDEF), where the energy 
spectrum was obtained from the first wall of the high neutron flux region of the GDT-FNS. The first 
wall neutron energy spectrum served as the external driver neutron source. 


2.4.3 Activity and specific activity of produced ”Mo 


To evaluate the efficiency of the *Mo production and usage rate of U, we defined the total 
activity of °Mo produced by the SBS in one day (24 h operation) as A [Ci/day], and the daily °Mo 
produced per unit mass of **°U as the specific activity (SA) [Ci/kg/day]. The amount of produced 
Mo is increased by U fission and reduced by its own decay. Thus, the number of Mo nuclides 
[N(t)] at time ¢ [s] changes according to equation (2) as follows: 


(2) 
where is the number of Mo nuclides at time ż [s], A is the decay constant of Mo, is the 
fission yield of ”Mo (0.061), = is the macroscopic fission cross-section of *°U [barns], and @ is the 
neutron flux tems ‘1. 
By defining and integrating equation (2), an activity equation (3) is obtained, which is 
consistent with the generalized activity equation[37]. 
(3) 
However, in an SBS, the activity equation can be modified to obtain equation (4) as follows: 
(4) 
In this case, the neutron flux (®) is replaced with , that is, the external neutron source intensity 
[n/s| from the high neutron flux region of the GDT-FNS. È is equal to , where [barns] is the average 
microscopic fission cross-section of a *°U atom. is the total number of *°U atoms in the SBS, 
which is equal to , where is the Avogadro constant (6.02x10”), M is the relative atomic mass of *°U 
equal to 235, and m is the mass of **°U[g] in the LEU solution 
In this study, was calculated by setting the tally 4 card and multiplier card FM4 in the input 
file. Therefore, by substituting the values into equation (4), the daily production of Mo in SBS 
can be evaluated. 


2.4.4 Neutron and gamma flux calculation 


The neutron and gamma flux can be calculated with a tally 4 [38] card by setting the different 
particle types [Neutron (N) or Photon (P)] using equation (5): 

(, (5) 
where is the neutron or gamma flux of the point detector [particles/cm’], (,,E,t) is the angular flux 
[particles/(cm*/sh/MeV/rad)], is the position vector [cm], E is the energy of the incident particle 
[MeV], is the direction vector, and is time [sh;lsh = 10 %s]. To evaluate the neutron flux 
distribution of the SBS, the heterogeneous coefficient Ky is introduced. This is defined as the ratio 
of the maximum value of the core thermal neutron flux to the average value. Ky is generally 
required to be less than 1.4 in a design [39]. 


2.4.5 Nuclear heat deposition 


The nuclear heat deposition was calculated using a tally 6 card (T6) combined with a 
superimposed mesh tally card (FMESH). The calculation result for T6 is the average energy 
deposition on the calculated cell, as shown in equation (6) [38]. 

(6) 
where is the total energy deposition in the cell [MeV/g], is the atom density (10%. atoms/cm’), m is 
the cell mass [g], is the microscopic total cross section [barns], and H(E) is the heating number 
[MeV/collision]. 


3 Results and Discussion 


3.1 Neutron Spectrum of the High Neutron Flux Region of GDT-FNS 


The neutron spectrum is a significant parameter of GD-FNS, which can affect the fission 
reaction efficiency. In this study, the neutron generation rate (Fig. 3) was used as the input data to 
calculate the neutron spectrum. The neutron spectrum of the plasma tube in the high-neutron- 
generation-rate region was obtained using the SuperMC program. The statistical error of the 
calculated neutron spectrum was ensured to be less than 1%. The spectral distributions are shown in 
Fig. 5. The results demonstrated that neutrons with an energy of approximately 2.5 MeV had the 
highest proportion, which is owing to the D-D reaction producing neutrons with an average energy 
of 2.5 MeV. In the subsequent neutronics design of the SBS, the calculated neutron spectrum is used 
as the external driver neutron source. 


Fig. 5 Neutron spectrum in the high neutron flux region 


3.2 Preliminary Design of the SBS 


The *Mo production by the SBS aims to achieve 50% of China’s projected °Mo demand by 
2025. China's current °°Mo heavily relies on imports, as demand continues to grow. China’s demand 
for medical ”Mo was estimated to be approximately 16000 6-day Ci in 2019 [40]. Considering an 
annual growth rate of 5% [5], the estimated °°Mo demand in 2025 will be approximately 21500 6- 
day Ci (i.e., 59 6-day per day). Determining these estimates requires the consideration of the decay 
losses during the separation and purification of the generated Mo. Approximately 80% of the 
originally produced ”Mo is lost during the 6-day period. However, approximately 10% [41] of the 
Mo cannot be extracted by chemical separation and purification processes. Therefore, the daily 
Mo demand is estimated to be 298 Ci in 2025 and 447 Ci in 2035. 

When determining the preliminary design of the SBS, the following constraints were 
considered: 1) To ensure the nuclear critical safety of the LEU solution SBS, k, must be less than 
0.98 [32-36]. However, considering the neutron source fluctuation and measurement uncertainty, as 
well as other uncertain factors, ks was limited to less than 0.97 to provide sufficient critical safety 
margins. 2) The U concentration was limited to range between 60-150 g/L. The upper limit, as 
indicated in the earlier sections, is owing to the saturation of UO2SO, at 150 g/L. 3) The SA should 
be as high as possible. 

The initial conditions for the calculation models were set as follows: the inner diameter of the 
plasma tube was 35 cm, the inner and outer radii of the LEU solution were between 37—74 cm, the 
thickness of the LEU solution container wall was 1 cm, the solid angle of the LEU solution was 
betweenz/4—n/3, and the thickness of the top and bottom reflectors was 5 cm. The back reflector 
was composed of a Be material with a thickness of 8 cm. There was a shielding thickness of 8 cm 
for materials composed of W, B, and polyethylene (PE) (Mass ratio 4:3:3). Different parameters of 
ks, A, and SA, were calculated by changing the U concentration and solid angle. The corresponding 
calculation results are listed in Table 3. 


Table 3 Preliminary calculation results of the different SBS design parameters (cases 1—6) 


235 
Case a d cen ae IL Bi ks A/Ci SA/Cikg' 
1 n/4 100 145.2 2.904 0.9465 92.5 31.9 
2 m/4 105 145.2 3.049 0.9671 130 42.6 
3 57/18 95 161.3 3.065 0.9410 91.7 29.9 
4 51/18 100 161.3 3.226 0.9612 125 38.7 
5 57/18 105 155.5 3.266 0.9681 156 47.8 
6 n/3 95 193.5 3.677 0.9576 143 38.9 


Note: A is the daily ”Mo production of one SBS; SA is the specific activity (daily °Mo production of 1 kg *°U). 


As shown in Table 4, the k, of case 5 satisfies the critical safety condition (less than 0.97), and 
the A is 156 Ci. This activity is close to the 50% medical *Mo demand for the projected Chinese 
market by 2025 (149 Ci). In addition, the SA of case 5 is the largest, indicating that the **°U use is 
also the most efficient compared to the other cases. Based on the results in Table 4, case 5 was 
selected as the preliminary scheme for the subsequent optimization design, including optimizing the 
U concentration, and the material types and sizes of the reflector and shielding. 


3.3 The Impact of Uranium Concentration on the SBS 


To study the influence of various U concentrations on the performance of the SBS *Mo 
production, the conditions of case 5 shown in Table 4 were selected for the detailed calculations 
with varying U concentrations. The volume of the LEU solution was set at 155.5 L, while the *°U 
enrichment was fixed at 19.75%. The distribution of the varying U concentrations, the 
corresponding changes in the neutron multiplication factor k,, and the daily Mo production 4 are 
shown in Fig. 6. 


Fig. 6 The impact of varying the U concentration (enrichment of 19.75% and solution volume of 
155.5 L) on the distribution of k, and A 

Based on the results shown in Fig. 6, there is a general increasing dependence of both k, and A 
on the increasing U concentration. In particular, the following can be deduced:1) As expected, the 
value of k, demonstrates a strong dependence on the U concentration because it sharply increases as 
the U concentration increases. However, note that at a U concentration of 110 g/L, the k, exceeds 
the design safety limit of 0.97. 2) Conversely, while the daily Mo production A increases gently 
with an increasing U concentration, there is a sharp increase when the U concentration is above 90 
g/L. 3) The most favorable condition in terms of the design nuclear critical safety is when the U 
concentration is just below 105 g/L, where k, is less than 0.97. Therefore, the U concentration was 
set to 105 g/L in the final design model. 


3.4 Optimal Design of Reflector 


The reflector is arranged on the outside of the LEU solution and can reflect the neutrons back 
into the U solution to minimize neutron leakage and improve the use of neutrons. The selection of 
the reflector material and size has an important influence on the efficiency of the ”Mo production 
and its associated k,. The commonly used reflectors are Be metal, BeO, graphite (GR), heavy water 


(D-H), zirconia (ZrO2), among others. In this section, case 5 is adopted for the calculation, and the 
thickness of the back reflector is set to 8 cm, while the other parameters remain unchanged. The 
values of ks were calculated by changing the reflector materials, the results of which are shown in 
Fig. 7. 


Fig. 7 Variations in k, using different reflector materials 


The k, values obtained using different reflector materials demonstrate a maximum value for 
the case of Be as a reflector, which indicates that Be presents the best reflector effect for the 
selected conditions. Therefore, Be was selected as the reflective material for the SBS design. The 
subsequent optimization (calculation) process for k, and A was performed by changing the thickness 
of the Be reflector. The results are shown in Fig. 8. 


Fig. 8 k,andA values varying with the thickness of the Be reflector 


The calculation results demonstrate that the values of k, and A increase as the thickness of the 
Be reflector increases. This trend emphasizes the beneficial effect of neutron use in SBS, with a 
corresponding benefit to the ”Mo production. In addition, the optimum thickness of the Be reflector 
is approximately 10 cm; beyond 10 cm, only a marginal increase in A is observed with a relatively 
significant increase in k.. Thus, to benefit from the capital cost of the reflector material, minimizing 
the geometric dimension of the entire SBS while improving the critical safety, a reflector thickness 
of 10 cm was plausible. 


3.5 Shielding Layer Optimization Design 


A shielding layer was used to reduce the neutrons and gamma radiation in the environment. 
Different shielding materials have varying shielding abilities for neutrons and gamma-rays. Eight 
types of materials were selected for the shielding study to reasonably select the material and 
thickness of the shielding layer; eight types of materials were chosen for the shielding study. The 
eight types of materials were Fe, Pb, W, Fe/B (weight ratio 1:1), W/B (weight ratio 1:1), Fe/PE 
(weight ratio 1:1), W/PE (weight ratio 1:1), and W/B/PE (weight ratio 4:3:3). The shielding 
performances of the different materials were evaluated with a thickness of 8 cm. The corresponding 
k, values and average neutron fluxes outside the shielding layer were calculated. The results and 
their distributions are shown in Fig. 9. 
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Fig. 9 Distributions of k, and the average neutron flux (outside the shielding layer) varying with 
different shielding materials 


The results demonstrate that the different shielding materials have little influence on the k, 
value, that is, the difference between the maximum (0.9695) and minimum values (0.9667) is only 
0.0028, as shown in Fig. 9. This indicates that W/B/PE had the best shielding effect for neutrons, 
which is expected because B has a very good thermal neutron absorption ability; PE is a good 


neutron moderator, while W, Fe, and Pb, have good gamma shielding. Overall, considering the 
shielding performance of the material against the neutrons and gamma rays and avoiding the toxic 


lead material, the composite material of W/B/PE was selected as the shielding material. 
The influence of the varying W/B/PE thicknesses on k, and the shielding performance were 


further studied. The calculation results are shown in Fig. 10, which demonstrate that k, presents no 
significant change as the shielding thickness increases because the shielding material contains B, 
which absorbs neutrons. However, the average neutron flux and gamma flux outside the shielding 
layer decrease with an increase in the thickness of the shielding material. 
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Fig. 10 The neutron flux and gamma flux outside the shielding layer and k, varying with the 
thickness of the W/B/PE shielding material 


According to the requirements of the shielding design, to limit the radiation impact from 
thermal neutron-activated products, the thermal neutron flux should be less than 1x10% cm”?s™, and 
[42, 43]. The results demonstrate that when the 
thickness of W/B/PE is 15 cm, the average neutron flux outside the shielding layer is 4.31<10* 


the gamma flux should be less than 4x10'° 


cms ! 


cms ', whereas the average gamma flux is 2.10x10’ cm’s' with k, = 0.9685. These values meet 
the requirements of shield design and nuclear critical safety margin. 


3.6 Neutron Flux Distribution of SBS 


The radial and axial neutron flux distributions of the SBS were calculated, the results of which 
are shown in Fig. 11 (a) and (b), respectively, demonstrating that the neutron flux in the U fission 
zone is of the order of 10" n/cm’:s (average of 3.73 x10" n/cm*. S; peak value of 4.94 x10" 
n/cm’:s). In addition, the neutron flux rapidly decreases after exiting the reflector layer (i.e., in the 
shielding layer). This sharp decrease confirms the effectiveness of the shielding materials indicated 
in the previous sections. In the axial direction, the neutron flux remains constant across the U 
solution, whereas the neutron flux at both ends of the reflector layer sharply decreases. The average 
neutron flux was calculated to be 3.88x10''n/cm’s and the maximum neutron flux was 
4.72x10''n/cm’s. As defined in Section 2.4.4, Ky can be calculated as 4.72<10''/3.88x10!! = 1.22, 
which meets the design requirements (Ky needs to be less than 1.4). This demonstrates that the 
radial and axial neutron flux distributions in the SBS are relatively uniform, which is beneficial for 
the efficient use of *°U and the safe operation of the system. 


Fig. 11 (a) Radial neutron flux distribution of SBS; (b) axial neutron flux distribution of SBS 


The neutron energy spectrum characteristics of the U solution layer, outer U container, reflector 
layer, and shielding layer were calculated as shown in Fig. 12. The results demonstrate that there are 
two thermal neutron peaks in the thermal neutron region (10 *-10°° MeV) and a fast neutron peak 
(approximately 2.5 MeV). The two thermal neutron peaks are owed to the H20 in the U solution, 
which served as a neutron moderator. This moderation causes many neutrons to be moderated into 
thermal neutrons. The fast neutron peak appears because the external neutrons driving the SBS are 
mainly the 2.5 MeV neutrons (D-D reaction neutrons). In addition, fast neutrons were produced by 
the fission of *°U. 


Fig. 12 Neutron energy spectrum of the different components in SBS 


3.7 Nuclear Heat Distribution of SBS 


To obtain the distribution of the nuclear heat in each component of the SBS, a tally card (T6) 
combined with FMESH was used to calculate the nuclear heat deposition. The visualization 
function of the SuperMC code was used to display the results, as shown in Fig. 13 (a) and (b). 


Heat deposition (W/cm}) Heat deposition (W/cm?) 


1. 55E-01 h ( ) 1. 55E-01 
1. 55E-03 1. 55E-03 
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1. 55E-07 1. 55E-07 
tenner Oe - 1. 55E-09 
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T 1. 55E-11 

1. 55E-13 
1. 55E-13 

1. 55E-15 
1. 55E-15 

1. 55E-17 
1. 55E-17 


(a) 


1. 55E-19 
1. 55E-19 


Fig. 13 (a) Nuclear heat distribution in the central radial section of SBS; (b) nuclear heat 
distribution in the central axial section of SBS 
The results demonstrate that the nuclear heat is mainly deposited in the LEU solution, with a 
maximum nuclear heat of 1.53x10 W/cm? at the central position. The minimum nuclear heat was 
determined to be 4.57x10 7° W/cm? at the edge position. The average nuclear heat was 7.01x107 
W/cm’, and the total nuclear heat was 10.9 kW. The average nuclear heat in the reflector layer was 
found to be 1.21<10* W/cm’, while the average nuclear heat in the shielding layer was 5.57x10° 
W/cm’. The nuclear heats of the reflector and shielding layers were approximately 2-3 orders of 
magnitude lower than that of the LEU solution. This is because the nuclear heat mainly arises from 
the fission energy generated by the fission of *°U, and the nuclear heat of the reflector and shielding 
layers is mainly from neutrons and gamma radiation energy deposition, which is significantly lower 
than the fission energy. COMSOL [44] was used to simulate the cooling system; according to the 
simulation results, the fuel solution may boil within 2 h if a cooling system is not added. The 
simulation results also demonstrate that a supplied inlet H2O coolant temperature of 22 °C and flow 
velocity of 1.0 m/s will be sufficient to maintain the fuel solution temperature below 90 °C. 
Based on the aforementioned analysis, the optimized design results for the °Mo production by 
SBS driven by GDT-FNS are listed in Table 4. 


Table 4 Final optimized design and parameters of Mo production SBS 


Length Thickness 


Zone (em) (en) Material 
Plasma at SBS 90 35 D 
Plasma tube 90 1 Fe 
LEU solution 90 36 UO.SO, aqueous solution, 
3U enrichment: 19.75%, U 
concentration 105g/L 
LEU solution container 90 1 Fe 
Reflector Back reflector 90 8 Be 
Top & bottom 54 5 
reflector 
Side reflector 90 8 
Shielding layer 90 15 W/B/PE 


Parameters Value 


k, 0.9685 
SA 48 Ci/daykg 
Yield of”Mo 157 Ci/day 


4. Conclusions 


In this study, an LEU SBS driven by GDT-FNS for ”Mo production was proposed. A 
neutronics analysis of the *Mo production system was conducted using the Monte Carlo method 
(SuperMC code). The neutronics analysis includes the calculation of the neutron spectrum in the 
region of the high neutron generation rate of the GDT-FNS. The analysis also covers the 
preliminary design and optimization assessments related to different U concentrations and their 
Mo production activities. Other analyses include the optimization design of the reflector and 
shielding layer, neutron flux, and the nuclear heat distribution of the SBS. 

In all optimization cases, the designed system must meet the safety requirements and the 
amount of Mo production necessary to meet 50% of China’s projected °Mo demand in 2025. The 
preliminary assessment, as shown in case 5 in Table 4, demonstrated the most favorable conditions, 
where the U concentration was 105 g/L for an LEU solution of 155.5 volume/L with a mass of 
3.266 kg and a subcritical multiplication factor of 0.9681. A further analysis was performed by 
placing the upper limit of the LEU solution volume at 155.5 and varying the mass of the LEU 
(19.75% enrichment) in the solution. The distribution was compared with the daily °°Mo production 
and its impact on the subcritical multiplication factor. Based on this analysis, the most favorable 
condition in terms of the designed nuclear critical safety was a U concentration of 105 g/L for a k, 
value near 0.97. Other analyses included shielding and reflector material selection and design 
optimization. Calculations demonstrated that Be (10 cm thick) and W/B/PE (15 cm thick) were 
suitable for serving as the reflector and shielding layers. The main optimized parameters are 
summarized as follows: 

1) The optimal value for the subcritical multiplication factor (ks) for the designed SBS was 
0.9685, while the average neutron flux and gamma flux outside the shielding layer were found to be 
4.31104 n/cm’s and 2.10x10’ n/cm’s, respectively. The distribution of the neutron flux and nuclear 
heating in the SBS were relatively uniform, as indicated by the Ky-value of 1.12, which further 
ensures the enhanced operational safety of the system. 

2) The SBS allows for a high **°U use, that is, 48 Ci ”Mo can be produced from 1 kg of *°U. 

3) A total of 157 Ci ”Mo can be produced by one SBS per day. Because the GDT-FNS is an 
axisymmetric structure and the solid angle of the SBS is only 57/18, multiple SBSs can be 
simultaneously arranged in the high neutron flux region of the GDT-FNS. According to the 
calculations, two and three of the SBSs designed in this study for the °Mo production can meet the 
demand of the Chinese market by 2025 and 2035, respectively. 

The SBS driven by the GDT-FNS ”Mo production system has the advantages of a high production efficiency, 
low nuclear waste, and low cost. Our study indicates that this system can be used as a potential facility for °Mo 


production. However, to make the system more feasible and practical, it is necessary to perform further detailed 
design studies, such as a U burnup analysis, ”Mo separation, and purification technique verification. 
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